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Abstract 

The  Single-Level  Breit-Wigner  method  of  microscopic  cross  section 
calculation  is  the  most  general  and  simplest  of  resolved  resonance  region 
formulas.  The  Single-Level  Breit-Wigner  formula  is  complicated  and  dependant 
on  multiple  parameters.  The  resonance  region  cross  section  formulas  require 
summations  and  multiple  trigonometric  function  calculations  to  arrive  at  cross 
section  values.  When  computer  codes  require  multiple  microscopic  cross  section 
values  inside  the  resonance  energy  region,  use  of  the  resonance  region  formula  to 
calculate  cross  sections  takes  considerable  computing  time.  In  an  effort  to  reduce 
computing  time  cross  section  values  are  often  approximated  by  piecewise  linear 
interpolation.  However,  the  use  of  piecewise  cubic  Hermit-e  polynomials  to 
approximate  resonance  region  cross  sections  provides  a  means  to  accomplish  the 
same  task  with  increased  accuracy  and  fewer  points  that  require  exact  cross 
section  calculation.  Once  the  derivative  of  the  cross  section  formula  is  derived 
and  the  interpolating  cubic  is  appropriately  shifted  and  scaled  for  numerical 
stability,  use  of  the  piecewise  cubic  Hermite  polynomial  is  easily  inserted  into 
computer  codes  that  depend  on  cross  section  calculation  in  the  resonance  energy 
region. 
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Efficient  and  Accurate  Computation  of  Elastic 


Cross  Sections  in  the  Single-Level  Breit-Wigner 
Resonance  Region 

I.  Introduction 

I.l:  Motivation 

The  neutron  cross  section  of  a  material  is  a  fundamental  piece  of 
information  for  many  important  applications.  For  instance,  radiation  health 
physics  and  reactor  physics  require  neutron  shielding  calculations  to  establish  the 
neutron  containment  capability  of  materials.  The  neutron  cross  section  is 
required  to  determine  a  containment  material’s  capability  to  scatter  or  slow  high 
energy  neutrons  and  also  absorb  neutrons  [2],  Reactor  physics  also  requires 
neutron  cross  sections  to  determine  the  criticality  conditions  of  a  reactor.  These 
criticality  calculations  require  neutron  cross  section  information  to  determine  a 
moderator’s  capability  to  scatter  and  slow  high  energy  neutrons  to  energies  more 
favorable  for  fission.  Additionally,  the  nuclear  cross  section  of  material  in  the 
reactor  core  in  various  geometries  is  required  to  determine  suitability  to  produce 
energy  by  fission  [3].  These  are  only  a  few  of  the  applications  that  depend  on 
neutron  cross  section  calculations  to  arrive  at  a  desired  conclusion.  Any  time 
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neutron  behavior  in  a  material  is  required,  neutron  cross  section  information  of 
the  material  is  required.  (See  Appendix  A  for  further  information  on  nuclear 
cross  sections). 

Calculation  of  multiple  microscopic  nuclear  cross  section  values  in  the 
resonance  region  of  a  material  is  computationally  expensive  in  terms  of  central 
processing  unit  (CPU)  time.  Computer  codes  that  calculate  cross  sections 
typically  approximate  cross  sections  in  the  resonance  region  by  piecewise  linear 
interpolation  of  the  cross  section  value  in  an  effort  to  make  efficient  use  of  CPU 
time.  Two  examples  of  such  codes  are:  The  Nuclear  Processing  Code  System 
(NJOY)  [4]  and  Piecewise-average  group  scattering  PAXK  [5].  NJOY  is  a 
production,  single  precision,  Fortran  77  code  that  is  sophisticated  but  difficult  to 
read  or  implement  into  research  codes.  PAXK,  developed  by  Gerts,  is  a  modern, 
double  precision,  Fortran  90/95  research  code  that  provides  accurate  cross  section 
calculations  and  is  well  modularized  and  straightforward  to  read  [5].  However, 
PAXk  fails  to  reach  the  accuracy  level  of  the  piecewise  linear  approximations 
produced  by  NJOY  with  a  reasonable  number  of  required  exact  calculations  that 
must  be  stored  to  produce  the  piecewise  approximation. 

1.2:  Goals  of  the  Research 

The  primary  goal  of  this  research  is  to  implement  a  piecewise 
approximation  to  cross  sections  in  the  resonance  region  with  the  following 
qualities: 
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1.  Accuracy  must  exceed  that  of  NJOY. 


2.  The  number  of  exact  calculations  stored  to  produce  the  piecewise 
approximation  must  be  fewer  than  required  by  NJOY. 

3.  The  implementation  must  be  relatively  simple  and  written  in  double 
precision,  modularized,  modern  Fortran  90/95  code. 

1.3:  Background 

The  resonance  energy  region  of  an  isotope  is  characterized  by  large 
changes  in  cross  section  values  with  corresponding  small  changes  in  incident 
neutron  energy.  Figure  1,  on  the  following  page,  shows  the  resonance  region  for 
Plutonium  238  (  238Pu)  at  zero  degrees  Kelvin  (K).  The  vertical  axis  depicts  the 
value  of  the  elastic  scattering  cross  section  in  barns  ( 10"28  meters2 )  and  the 
horizontal  axis  depicts  incident  neutron  energy  in  electron  volts  (eV).  Figure  2  is 
a  more  isolated  view  of  the  2.855  eV  resonance  peak.  The  vertical  axis  of  both 
figures  is  presented  on  a  logarithmic  scale.  Figure  2  shows  the  difficulty  in 
creating  a  piecewise  linear  approximation  to  cross  sections  in  the  resonance 
region  by  the  displaying  the  obvious  nonlinear  behavior  of  the  resonance  peaks. 
The  cross  section  calculations  used  to  construct  the  plots  in  Figure  1  and  Figure 
2,  were  conducted  in  energy  increments  of  1  x  10"5  eV. 
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The  plots  in  Figure  1  and  Figure  2  were  generated  by  PAXK  calculation  of 

the  Single-Level  Breit-Wigner  (SLBW)  formula  across  the  resonance  region  using 

parameters  found  in  File  2  of  2,i8Pu’s  associated  Evaluated  Nuclear  Data  File 

(ENDF)  (see  Appendix  B  for  more  detailed  information  on  ENDF).  The  ENDF 

provides  parameterization  and  recommended  calculation  methods  to  construct 

the  resonance  region  cross  sections.  This  ENDF  information  produces  cross 

section  values  at  zero  K  (NJOY  and  PAXK  perform  Doppler  broadening 

calculations  to  construct  cross  sections  at  other  temperatures). 

The  SLBW  formula  and  first  expansion  are  presented  below  to  indicate 

why  an  interpolating  polynomial  of  the  form  c0x  +  q  or  even 
o  o 

C3X  +  qaT  +  qx  +  c0  is  computationally  advantageous.  The  SLBW  elastic 
scattering  formula  is, 

nls— 1 

<7„(£)  =  £<7,u(£).  (1) 

1=  0 


The  first  expansion  of  Equation  (1)  is, 


<7u  (E)  =  (2(  +  1)777777  sin2  (<Pi  ( E ))  + 


E  (E) 


7 r 


NR, 


E  sj  E 


U[(E-E'r{E))  +YAT2r{E) 

[K,r  (E) COS (2^  (E))  -2Tnr  ( E)Txr  sin2  (0;  (E))  + 
2[E  —  E'r  (E))  rn>r  sin  (20;  (E))) 


(2) 


(See  Appendix  C  for  explanation  of  terms  and  full  expansion  of  SLBW  formula.) 
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As  stated  earlier,  NJOY  uses  undocumented  sophisticated  techniques  to 


reconstruct  the  resonance  region  in  a  piecewise  linear  fashion.  Documentation 
indicates  that  the  initial  points  stored  for  construction  of  piecewise  linear 
interpolation  are  energies  of  resonance  peaks  or  “half-height”  energies  of 
resonance  [6].  The  energies  of  the  resonance  peaks  are  found  in  the  ENDF 
parameter,  resonance  energy  ( Er )  [7].  “Half-height”  resonance  energies  are  only 
mentioned  once  in  the  NJOY  documentation  and  are  not  defined.  Convergence 
criteria  and  methods  for  selection  of  the  final  number  and  location  of  exact 
calculations  that  require  storage  for  the  piecewise  linear  approximation  of  cross 
sections  in  the  resonance  region  is  not  stated  explicitly.  The  criteria  and 
methods  are  only  characterized  as  more  “complicated”  than  those  used  for 
energies  in  the  smooth  regions  [6].  Whatever  techniques  are  used  by  NJOY  to 
establish  the  piecewise  linear  approximation  are,  they  produce  more  accurate 
results  than  those  of  PAXK. 

The  PAXk  method  of  producing  a  piecewise  linear  approximation  to 
resonance  region  cross  sections  is  much  more  straight-forward.  The  Er 
parameters  are  input  and  corresponding  exact  cross  section  calculations  are 
made.  The  same  procedure  is  completed  for  the  endpoints  of  the  resonance 
region.  These  endpoints  and  Er  calculations,  are  then  stored  as  the  initial  set  of 
required  calculations  for  piecewise  interpolation.  PAXK  then  applies  an  adaptive 
midpoint  method  to  generate  the  required  energy  mesh  for  piecewise  linear 
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approximation.  Starting  with  the  midpoint  of  the  left  endpoint  and  first  Er 
value  desired  tolerance  checks  are  made  between  the  midpoint  approximation  and 
the  exact  cross  section  value  at  the  midpoint.  If  tolerance  is  met,  no  further 
energies  of  exact  calculation  are  stored.  If  tolerance  is  not  met,  the  quarter  point 
(that  is,  the  point  half  way  between  the  original  midpoint  and  the  left  endpoint) 
closest  to  the  first  end  point  is  checked  for  tolerance.  If  tolerance  is  met,  the  mid 
point  is  then  added  as  a  required  point  for  piecewise  interpolation.  If  tolerance  is 
not  met,  the  process  repeats  moving  closer  to  the  left  end  point.  This  procedure 
continues  until  the  midpoint  calculated  using  the  right  end  of  the  resonance 
region  is  accepted  as  a  required  point  for  piecewise  interpolation.  (See  Appendix 
E  for  an  algorithm  of  this  process). 

Although  the  PAXK  techniques  are  easy  to  understand,  the  results  are  not 
as  accurate  or  efficient  as  those  of  NJOY.  Comparisons  show  that  PAXK  requires 
an  increase  in  required  linear  approximations  by  a  factor  of  two,  to  approach  the 
accuracy  of  NJOY  piecewise  interpolation. 

1.4:  Problem  Statement 

Efficiently  and  accurately  reconstructing  resonance  region  cross  sections  is 
a  problem  for  modern  research  computing  codes,  such  as  PAXK.  The  use  of  the 
piecewise  cubic  Hermite  approximation  as  opposed  to  linear  approximation 
presents  a  viable  solution. 
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1.5:  Scope 


For  reasons  of  generality  and  simplicity,  the  Single-Level  Breit-Wigner 
(SLBW)  resonance  region  formula  is  chosen  as  the  method  for  resonance  region 
cross  section  calculation  for  this  research.  Although  more  sophisticated  formulas 
exist,  such  as  Multi-Level  Breit-Wigner  (MLBW),  Reich-Moore,  Adler- Adler  and 
Hybrid  R- function,  the  use  of  SLBW  is  widespread  in  engineering  codes  [7]. 

ENDF  documentation  provides  methods  to  convert  parameters  for  the  use  of 
SLBW  in  place  of  all  other  formulas  for  cross-section  calculation  in  resolved 
resonance  regions.  In  the  case  of  MLBW  and  Hybrid  R-function  the  parameters 
are  an  exact  substitute  for  use  in  SLBW  [7]. 

238 Pu  is  chosen  as  the  target  element  for  cross  section  calculations.  It  is 
one  of  the  few,  and  possibly  the  only,  elements  in  the  ENDF  library  that  actually 
has  a  resonance  region  where  the  SLBW  method  of  calculation  is  specified  for  use 
by  ENDF.  Also,  the  dominant  neutron  interaction  in  the  238 Pu  SLBW  resonance 
range  is  elastic  scattering. 

The  SLBW  formula  considers  three  neutron/nucleus  interactions:  Elastic 
scattering,  radiative  capture  and  fission.  To  reduce  the  difficulty  in  debugging 
the  computational  methods,  only  the  elastic  scattering  neutron  cross  section  (crel) 
is  calculated.  The  elastic  cross  section  formula  (see  equations  (1)  and  (2))  is  the 
most  mathematically  complicated  neutron  interaction  formula  in  the  SLBW 
method[7].  Because  of  the  complexity  of  the  elastic  scattering  cross  section 


formula,  any  interpolation  technique  that  is  successful  for  elastic  scattering  is 
expected  to  be  easily  expanded  to  other  interactions. 

1.6:  Approach 

Research  of  NJOY  documentation  provides  insight  into  the 
fundamental  requirements  of  approximating  the  resonance  region  [6].  These 
fundamental  requirements  are:  Calculation  of  an  initial  energy  mesh,  adaptive 
calculation  of  the  final  energy  mesh  and  piecewise  interpolation  between  mesh 
points. 

The  proposed  method  of  initial  energy  mesh  calculation  is  to  include  in  the 
initial  energy  mesh  the  energy  values  of  the  endpoints  of  the  resonance  region, 
and  the  energies  of  minimum  and  maximum  cross  section  values  at  each  of  the 
resonance  peaks.  The  proposed  method  of  adaptive  final  mesh  calculation  is  a 
slightly  modified  version  of  the  adaptive  midpoint  method  used  by  PAXK.  The 
proposed  method  of  interpolation  between  energy  mesh  points  is  piecewise  cubic 
Hermite  interpolation.  These  fundamental  requirements  and  a  theoretical 
analysis  of  the  suitability  of  the  proposed  techniques  are  presented  in  more  detail 
in  Chapter  Two. 

The  specific  computational  techniques  applied  to  meet  the  fundamental 
requirements  of  resonance  region  cross  section  approximation  are  presented  in 
Chapter  Three.  Derivations  required  for  implementation  of  the  computational 
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techniques  of  Chapter  Three  are  available  in  Appendix  C  and  Appendix  D. 
Validation  of  the  applied  computational  techniques  is  presented  in  Chapter  Four. 

A  statistical  error  analysis  detailing  the  increased  accuracy  of  the  applied 
computational  techniques  over  that  of  NJOY  is  found  in  Chapter  Five.  Chapter 
Five  also  includes  a  comparison  of  CPU  efficiency  between  linear  approximation, 
cubic  Hermite  approximation  and  SLBW  calculation  of  resonance  region  cross 
sections.  This  comparison  is  used  to  motivate  the  computational  advantage  of 
cross  section  approximation  when  compared  to  SLBW  calculations.  The  CPU 
timing  comparison  also  quantifies  the  trivial  computational  cost  of  using  cubic 
Hermite  interpolation  as  opposed  to  linear  interpolation. 
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II.  Theory 


To  produce  a  suitable  piecewise  approximation  to  the  resonance  cross 
sections,  three  fundamental  properties  are  required.  First,  the  function  must  be 
initially  characterized  to  avoid  inaccurate  approximations.  This  first  process 
includes  selectively  choosing  points  for  exact  calculation  and  storage  that  are 
necessary  to  avoid  anticipated  inaccurate  interpolation.  Second,  an  adaptive 
method  for  storing  only  enough  points  and  exact  calculations  required  to  meet 
desired  accuracy.  Third,  an  accurate  method  of  interpolation  is  required  to 
approximate  the  cross  sections  between  stored  points.  The  following  sections 
discuss  these  fundamental  properties  as  applied  to  SLBW  resonance  region  cross 
section  approximation. 

II.  1:  Initial  Characterization  of  the  SLBW  Resonance  Region  Formula 
As  stated  in  chapter  one,  the  resonance  region  is  characterized  by 
resonance  peaks.  In  the  case  of  the 238 Pu  SLBW  resonance  region,  there  are 
fourteen  peaks  in  the  region.  Figure  3  depicts  the  SLBW  calculation  and  the 
corresponding  NJOY  linear  interpolation.  Figure  4  illustrates  the  improvement 
in  accuracy  when  the  maximum  value  of  the  resonance  peak  is  included  in  the 
linear  approximation.  Similar  accuracy  increases  occur  when  corresponding 
minimum  values  are  included  in  the  approximation.  For  this  reason,  the 
proposed  method  stores  the  local  maximum  and  minimum  energies  of  resonance 
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peaks  to  initially  characterize  the  SLBW  resonance  region.  These  stored  points 


serve  as  a  starting  point  for  the  adaptive  methods  of  the  next  section. 


. 5615  18.5618  18.562  18.5623  18.5625  18.5628  18.563 

energy  (el?) 


Figure  3.  NJOY  approximation  to  the  Pu  18.56  eV  resonance  peak 
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Figure  4.  Linear  approximation  to  the  238Pu  18.56  eV  resonance  peak  with  local 
maximum  included 


12 


II.2:  Adaptive  Calculation  of  the  Energy  Mesh 


The  energy  mesh  is  the  collection  of  energy  values  and  associated  exact 
cross  section  calculations  that  require  storage  to  produce  piecewise  interpolation. 
The  use  of  an  adaptive  mesh  generator  greatly  reduces  the  number  of  points  in 
the  energy  mesh  by  only  keeping  those  points  which  are  required  to  meet  desired 
accuracy.  Referring  to  Figure  5,  points  A,  B  and  C  are  separated  in  energy  by  .1 
eV.  An  acceptable  interpolation  between  points  A  and  B  could  conceivably  be 
achieved  by  a  single  line.  This  is  most  definitely  not  the  case  for  interpolation 
between  B  and  C  if  similar  accuracy  is  desired.  More  points  are  required  to 
accurately  interpolate  between  points  B  and  C. 


Figure  5  .  236Pu,  2.855  eV  resonance  peak 
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Uniformly  dividing  the  energy  range  of  Figure  5  (2.7  eV  to  3  eV)  until  an 
acceptable  piecewise  approximation  was  achieved  would  result  in  unneeded  points 
between  points  A  and  B  of  Figure  5.  For  this  reason  an  adaptive  technique,  such 
as  the  bisection  tolerance  technique  used  in  PAXK,  is  required  to  minimize  stored 
points  in  the  energy  mesh. 

II. 3:  Piecewise  Cubic  Hermite  Polynomial  Approximation 

As  previously  mentioned,  both  NJOY  and  PAXK  use  a  piecewise  linear 
technique  to  approximate  the  resonance  cross  sections.  Piecewise  linear 
approximation  has  a  disadvantage  in  that  there  is  no  differentiability  at  the 
endpoints  of  the  subintervals.  This  translates  to  a  representation  that  is  not 
“smooth”  in  a  geometrical  sense  [8].  The  cubic  Hermite  polynomial  has  a 
continuous  derivative  across  the  energy  mesh  resulting  in  a  “smooth” 
representation  [8].  This  cubic  Hermite  approximation  meets  the  physical 
conditions  of  the  SLBW  resonance  formula  much  better  than  that  of  the  linear 
representation.  The  grayscale  line  of  Figure  6,  on  the  following  page  is  a  dense 
list  plot  of  exact  SLBW  calculations.  The  dashed  line  is  the  NJOY  linear 
approximation  and  the  black  line  overlaying  the  grayscale  is  the  cubic  Hermite  fit 
using  the  same  endpoints  as  the  NJOY  approximation.  More  detailed 
information  on  calculation  of  the  cubic  Hermite  fit  is  presented  in  the  next 
chapter.  It  is  clear  from  Figure  6  that  the  potential  for  increased  accuracy  is 
significant  when  using  a  cubic  interpolation  instead  of  a  linear  interpolation.  For 


14 


now,  it  is  assumed  that  the  benefits  of  increased  accuracy  and  fewer  required 


points  in  the  energy  mesh  outweigh  the  computational  costs  of  interpolating 
using  a  cubic  calculation,  C3X’  C2X’  +  °ix  +  co . 


en«rgy(«V) 


SLEW 

N  JOY 
Cubic 


Figure  6.  Hermite  Cubic  interpolation  of  the  238 Pu  18.56  eV  resonance  peak 
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Ill  Methodology 


Development  of  the  following  techniques  for  SLBW  resonance  cross  section 
approximation  was  accomplished  using  the  ENDF  processing  and  SLBW 
resonance  region  cross  section  calculation  methods  of  PAXK.  All  code  used 
during  the  research  was  written  in  Fortran  90/95  and  compatibility  with  original 
PAXk  code  was  maintained.  However,  each  of  the  following  three  sections 
describes  a  single  subroutine  or  module  that  is  fully  independent  of  the  other  two. 
The  methods  described  are  general  and  any  of  the  three  could  be  implemented  or 
improved  upon  separately.  The  only  unique  requirements  are:  calculation  of  the 
first  derivative  of  the  SLBW  formula  for  use  of  either  the  initial  energy  mesh 
calculator  or  the  piecewise  cubic  Hermite  approximation;  and,  an  initial  energy 
mesh  of  at  least  two  points  is  required  for  use  of  the  adaptive  energy  mesh 
calculator.  This  independence  among  computational  modules  is  important  to 
allow  further  improvements  in  the  presented  techniques  to  be  made  on  any  of  one 
of  the  techniques  with  minimal  impact  on  the  other  two.  It  also  allows  authors 
of  research  code  to  pick  any  of  the  three  methods  for  implementation,  at  his  or 
her  discretion,  without  being  forced  to  use  the  other  two. 

III.l:  Calculation  of  the  Initial  Energy  Mesh 

The  Er  parameters  in  ENDF  are  used  in  calculation  of  SLBW  resonance 
region  cross  sections.  The  positive  Er  values  correspond  to  the  energy  of  a 
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particular  resonance  peak  in  the  region.  Inspection  of  Er  values  provided  in 
ENDF  indicates  that  these  values  are  in  the  vicinity  of  resonance  peaks  but  do 
not  provide  the  exact  energy  location  of  the  peak.  A  more  precise  energy 
location  is  determined  by  use  of  the  SLBW  first  derivative.  (The  first  derivative 
of  the  SLBW  formula  is  a  significant  contribution  to  both  the  calculation  of  the 
initial  energy  mesh  and  the  piecewise  cubic  Hermite  approximation  of  section 
III. 3.  The  derivation  of  the  SLBW  first  derivative  is  provided  in  Appendix  C.) 

The  calculation  of  the  maximum  of  a  resonance  peak  begins  by  using  Er 
parameters  as  starting  points,  then  subtracting  and  adding  a  AE  value  of 
increasing  magnitude  until  the  signs  of  the  derivative  of  the  Er  —  AE  and  the 
derivative  of  Er  +  AE  are  opposite.  Once  this  condition  is  met,  a  bisection 
search  root  solver  approximates  the  energy  location  where  the  first  derivative  of 
the  SLBW  formula  is  equal  to  zero  and  thus  the  maximum  of  the  resonance  peak. 
The  approximation  is  accomplished  by  setting  the  convergence  criteria  of  the  root 
solver  to  an  absolute  difference  of  1CT14  between  the  current  and  the  previous 
approximated  energy  value.  The  calculation  of  the  corresponding  minimum  value 
of  a  resonance  peak  is  accomplished  similarly  with  the  following  exception.  The 
starting  point  isEr  —  AE  .  The  ENDF  defined  endpoints  of  the  resonance  region 
and  these  values  of  resonance  peak  maximum  and  minimum  values  are  then 
placed  into  the  initial  energy  mesh  (See  Appendix  E  for  algorithm). 
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SLBW  cross  section  calculations  of  approximated  resonance  peak  maxima 


and  minima  were  compared  to  those  of  energy  values  on  either  side.  Table  1, 
contains  values  of  the  resonance  peak  maxima,  minima  and  corresponding  Er 
parameter  for  the  238 Pu  SLBW  resonance  region. 


MINIMA  (EV) 

MAXIMA  (EV) 

Er  (EV) 

2.8299068750 

2.8694042930 

2.855 

9.9467878721 

9.9872619861 

9.975 

18.359319616 

18.562246850 

18.560 

59.744456589 

59.806950102 

59.800 

70.020724600 

70.106396134 

70.100 

82.484396635 

83.001647414 

83.000 

109.87398263 

110.00471228 

110.000 

113.56253961 

113.80270119 

113.800 

117.90064368 

118.60165448 

118.600 

121.93378776 

122.50207828 

122.500 

150.72762650 

151.00425399 

151.000 

169.98246548 

171.20204065 

171.200 

181.91672822 

182.50295768 

182.500 

191.14852354 

192.00326405 

192.00 

Table  1:  238  Pu  SLBW  resonance  region  peak  max,  min  and  ER  values 


The  approximated  values  of  the  local  maxima  and  minima  energy 
locations  proved  accurate  to  a  minimum  of  ten  digits  of  precision.  Comparisons 
made  beyond  the  tenth  digit  of  precision  returned  equal  SLBW  calculated  cross 
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sections  for  low  energy  resonance  peaks.  Since  the  SLBW  cross  section  values 
were  equal  for  energy  input  that  differed  beyond  the  tenth  digit  confirmation  of 
the  energy  location  of  the  minimum  and  maximum  values  of  resonance  peaks 
cannot  be  confirmed  beyond  the  tenth  digit. 

III. 2:  Implementation  of  the  Adaptive  Energy  Mesh  Calculator 

The  adaptive  energy  mesh  calculator  is  an  improved  version  of  the 
bisection  method  used  by  PAXK.  The  adaptive  energy  mesh  calculator  takes  the 
initial  energy  mesh  as  input  and  calculates  the  SLBW  cross  section  and  the 
piecewise  approximation  of  the  midpoint  between  the  first  and  second  value  in 
the  initial  energy  mesh.  A  user-specified  tolerance  check  is  then  made  at  the 
midpoint  of  both  absolute  and  relative  error,  where  absolute  error  equals 

|SLBW  -  Approximation!  and  relative  error  equals  SLBVV  -  Approximation  Jf 

SLBW 

the  tolerance  is  met,  no  further  points  are  added.  If  the  specified  tolerances  are 
not  met,  the  second  value  in  the  initial  mesh  is  temporarily  replaced  by  the  mid 
point  and  the  process  repeats  until  a  temporary  point  meets  tolerance  and  is 
accepted  into  the  energy  mesh.  The  entire  process  then  begins  again  with  the 
newly  added  mesh  point  and  the  second  initial  mesh  point.  Once  the  midpoint 
between  the  most  recently  added  mesh  point  and  the  second  initial  mesh  meets 
specified  tolerance,  the  process  begins  again  with  the  second  and  third  initial 
mesh  point.  The  process  continues  across  the  resonance  region  until  the 
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midpoint  between  the  most  recently  added  mesh  point  and  the  resonance  region 
endpoint  meets  tolerance.  (See  Appendix  E  for  algorithm). 

It  is  important  to  note  that  the  specified  tolerance  is  not  the  guaranteed 
accuracy  of  all  possible  approximations.  Results  presented  in  Chapter  4  show 
that  in  most  cases  accuracy  is  better  than  that  specified,  but  non-typical  error  is 
probable  for  some  sub-regions  of  the  resonance  region. 

III. 3:  Implementation  of  the  Piecewise  Cubic  Hermite  Polynomial 

The  piecewise  cubic  Hermite  Polynomial  is  implemented  in  much  the  same 
manner  as  piecewise  linear  approximation.  Both  methods  use  sequential  sets  of 
two  points  in  the  energy  mesh  and  corresponding  cross  section  values  to 
interpolate  values  in  between  the  two  points.  The  piecewise  cubic  Hermite 
polynomial  requires  the  additional  value  of  the  SLBW  first  derivative  at  the  two 
mesh  points.  The  value  of  the  cross  section  and  the  first  derivative  at  each  of 
these  two  mesh  points  combine  to  form  the  required  system  of  four  equations  to 
solve  for  the  four  unknown  coefficients  of  the  cubic  Hermite  polynomial, 
c3x’3  +  c2x2  +  cxx  +  c0  [8]. 

For  the  purposes  of  approximating  the  SLBW  resonance  region  of  Pu238, 
the  Hermite  cubic  was  scaled  and  shifted  by  a  change  of  variable.  This  was  done 
in  order  to  avoid  overflow  and  catastrophic  cancellation  in  the  coefficients. 
Accuracy  loss  on  the  energy  value  due  to  the  change  of  variable  proved 
insignificant  when  compared  to  linear  approximations.  Results  of  this 
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comparison  are  presented  in  Chapter  4.  The  piecewise  cubic  Hermite  coefficients, 


derivation  of  the  shifted  and  scaled  coefficients,  and  a  discussion  on  possible 
accuracy  loss  is  presented  in  Appendix  D. 
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IV  Validation  of  Computational  Methods 


IV.  1:  Validation  of  the  Initial  Energy  Mesh  Calculator 

The  module  responsible  for  calculation  of  the  initial  energy  mesh  relies  on 
an  accurate  calculation  of  the  first  derivative  of  the  SLBW  formula.  The  derived 
first  derivative,  found  in  Appendix  C,  was  initially  verified  by  calculating  a  dense 
series  of  energies  and  SLBW  cross  section  values  along  with  corresponding  SLBW 
first  derivative  calculations.  The  energies  used  for  these  calculations  were  sorted 
from  lowest  to  highest  and  were  separated  by  lxlO"5  eV  across  the  entire  238Pu 
SLBW  resonance  region.  The  SLBW  cross  section  calculations  at  the  first  and 
third  points  of  energy  in  a  series  of  three  energy  points  were  used  to  generate  a 
linear  slope  between  the  calculated  cross  section  values.  This  linear  slope  was 
then  compared  to  the  calculated  derivative  at  the  second  point  in  the  three  point 
series.  Relative  error  between  the  linear  slope  and  the  calculated  first  derivative 
was  less  than  0.002  across  the  entire  resonance  region.  Further  validation  of  the 
first  derivative  calculation  and  the  initial  energy  mesh  calculator  is  provided  in 
the  accuracy  of  the  local  minimum  and  maximum  energy  locations  calculated  by 
the  methods  described  in  chapter  three  with  results  provided  in  Table  1.  These 
energy  values  of  the  local  maximum  and  minimum  cross  sections  were  verified  by 
comparing  cross  section  calculations  on  either  side  of  the  calculated  energy 
locations  of  the  local  maximum  and  minimum  cross  sections.  As  stated  in 
chapter  three  these  values  proved  accurate  to  the  tenth  digit  of  precision.  This  is 
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not  to  say  that  the  calculated  energy  location  of  local  maximum  and  minimum 
cross  section  values  is  not  accurate  beyond  the  tenth  digit.  It  is  just  difficult  to 
verify  accuracy  beyond  the  tenth  digit  due  to  the  fact  that  a  change  in  energy 
beyond  the  tenth  digit  requires  more  than  double  precision  to  display  a  change  in 
cross  section  value. 

IV.  2:  Validation  of  the  Adaptive  Final  Energy  Mesh  Calculator 
The  concept-  of  midpoint  adaptive  final  energy  mesh  calculator  was 
previously  implemented  in  PAXK.  For  this  reason,  only  basic  debugging  was 
conducted  to  ensure  that  initial  energy  mesh  points  were  input  properly  and 
proper  values  were  captured  in  the  final  energy  mesh  without  duplicate  entries. 
Validation  was  limited  to  inspection  of  the  final  energy  mesh  to  ensure  expected 
regions  of  low  density  mesh  points  as  well  as  regions  of  expected  high  density 
mesh  points  were  present.  Figure  7,  shows  a  list  plot  of  the  final  energy  mesh 
calculated  using  the  methods  described  in  chapter  three.  The  vertical  axis  of 
Figure  7  is  logarithmic. 
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Figure  7.  238  Pu  SLBW  energy  mesh 

IV.  3:  Validation  of  Piecewise  Cubic  Hermite  Interpolation 

Piecewise  cubic  Hermite  polynomial  approximation  requires  an  accurate 
first  derivative  of  the  SLBW  formula  to  provide  accurate  interpolation.  The 
accuracy  of  the  SLBW  first  derivative  was  determined  during  validation  of  the 
initial  energy  mesh  calculator.  Further  validation  of  the  cubic  Hermite 
interpolating  polynomial  calculations  are  provided  by  the  results  in  the  following 
chapter. 
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V  Results 


To  determine  validity  of  the  resonance  region  cross  section  approximation 
techniques  presented  in  Chapter  3  (hereafter  referred  to  as  the  test  method),  a 
comparison  was  made  to  the  approximation  produced  by  the  production  ENDF 
processing  code  NJOY.  The  comparison  considered  total  number  of  energy  mesh 
points,  and  associated  data  storage  per  point,  required  to  accomplish  the 
resonance  region  approximation  of  238  Pu  as  well  as  a  statistical  comparison  of 
error.  Additionally,  code  profiling  was  conducted  in  order  to  establish  the 
computational  cost  of  using  piecewise  cubic  Hermite  approximations  as  opposed 
to  linear  approximations. 

PAXk  resonance  region  cross  section  approximation  is  not  presented  for 
detailed  comparison.  The  bisection  adaptive  mesh  calculator  when  using  linear 
interpolation  requires  an  increase  in  mesh  points  by  a  factor  of  two  in  order  to 
come  within  two  orders  of  magnitude  of  the  mean  relative  error  of  NJOY.  For 
this  reason  it  is  excluded  from  the  following  results  in  order  highlight  the  more 
meaningful  comparison  of  the  test  method  with  NJOY  resonance  region  cross 
section  approximation. 

IV.  1:  Analysis  of  Data  Storage  to  Approximate  the  Cross  Section 

The  NJOY  energy  mesh  for  the  Pu238  SLBW  resonance  region  contains 
6360  points,  containing  two  values  each  (Fand  crd(E))  resulting  in  12,720  stored 
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values  [6].  The  energy  mesh  using  piecewise  cubic  Hermit-e  approximation 
required  1380  energy  mesh  points,  with  a  specified  tolerance  of  0.01  barns  of 
absolute  error  and  a  relative  error  of  0.0001.  Each  mesh  point  contains  three 
values  ( E  ,  a(E )  and  a7  (A) )  resulting  in  4140  stored  values.  The  test  method 
reduced  the  number  of  energy  mesh  points  by  a  factor  of  four  and  reduced  the 
data  storage  requirement  of  the  energy  mesh  by  a  factor  of  three.  The  economy 
and  accuracy  of  energy  mesh  points  of  the  two  methods  is  illustrated  in  the 
following  figures. 

Figure  8,  shows  the  SLBW  calculation  in  gray  along  with  the  piecewise 
cubic  Hermite  approximation  and  the  NJOY  approximation.  The  cubic  Hermit-e 
approximation  is  presented  using  calculations  of  the  shifted  and  scaled 
polynomial  discussed  in  Chapter  3  and  Appendix  D.  Both  the  NJOY  and  test 
approximation  used  three  energy  mesh  points  to  approximate  the  cross  section  in 
this  region.  It  is  apparent  from  the  figure  that  NJOY  would  require  many  more 
points  to  approximate  the  cross  section  as  accurately  as  the  test  method. 
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Figure  8  .  238Pu  182.5  eV  resonance  peak  approximation 


Figure  9,  shows  the  NJOY  approximation  to  the  cross  section  in  the 
vicinity  of  the  local  minimum  value  of  the  Pu238  182.5  eV  resonance  peak. 

NJOY  used  seven  energy  mesh  points  to  approximate  this  region.  It  is  seen  from 
Figure  9  that  NJOY  made  a  good  approximation  of  the  cross  section  in  this  case. 
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Figure  9  .  238Pu  182.5  eV  local  minimum  NJOY  approximation 
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Figure  10,  adds  the  test  methods  approximation  to  the  plot.  The  test 


method  used  only  three  energy  mesh  points  to  approximate  cross  section  in  this 
region  to  a  better  accuracy  than  that  of  the  NJOY  approximation. 
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Figure  10  .  238 Pu  182.5  eV  local  minimum  approximation 


IV.  2:  Error  Analysis  of  The  Cross  Section  Approximations 


The  statistical  error  analysis  consists  of  SLBW  calculations  of  100,000 
randomly  generated  points  of  energy.  The  calculations  are  compared  with  the 
NJOY  approximation  and  the  test  method  approximation  for  each  of  the  100,000 
points  of  energy.  The  error  analysis  includes  both  the  absolute  and  relative  error, 
where,  absolute  error  equals  |SLBW  -  Approximation)  and  relative  error 
|  SLBW  -  Approximation  | 


equals 


SLBW 


.  The  energy  value  of  the  NJOY  energy  mesh 


was  produced  directly  by  NJOY  in  single  precision.  The  corresponding  cross 
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section  value  was  calculated  by  PAXK  in  double  precision.  The  double  precision 


calculation  reduced  the  error  of  the  NJOY  approximation  and  was  used  to  make 
comparisons  with  the  test  method  approximations  as  equitable  as  possible. 

The  mean  error  of  both  the  NJOY  approximation  and  the  test  method 
approximation  is  less  than  the  corresponding  standard  deviations.  Because  of 
this  inconclusive  combination  of  mean  and  standard  deviation,  percentile  plots  of 
error  are  presented  in  following  figures.  Percentile  plots  depict  the  error  on  the 
vertical  axis  and  the  percentage  of  total  data  points  along  the  horizontal  axis. 
The  plotted  value  indicates  what  percentage  of  data  points  fall  below  the 
indicated  error  value.  Percentile  plots  give  an  overall  characterization  of  the 
typical  and  non-typical  error.  The  proceeding  figures  indicate  error  that  is 
influenced  by  few  relatively  large  non-typical  results  in  both  the  case  of  absolute 
and  relative  error.  For  ease  of  viewing  only  the  greatest  one  percent  of  the  data 
points  is  depicted  on  the  percentile  plots.  The  vertical  axis  includes  the  data 
point  of  the  greatest  error. 
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Figure  11.  Percentile  plot  of  NJOY  absolute  error,  99-100% 


Figure  12.  Percentile  plot  of  Test  Method  absolute  error  99-100% 
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Figure  13.  Percentile  plot  of  NJOY  relative  error  99-100% 


0.006  F-r 


0.005 


0.004 


0.003 


0.002 


0.001 


99  99. Z  99.4  99.6  99. S 

P  e  r  c  tut 


Figure  14.  Percentile  plot  of  Test  Method  relative  error  99-100% 


The  cause  of  these  non-typical  error  results  in  the  case  of  NJOY  is  not 


known.  However,  for  the  test  case  it  is  attributed  to  the  tolerance  criteria  of  the 
adaptive  mesh  calculator.  The  typical  energy  location  of  non-typical  absolute 
error  is  in  the  vicinity  of  the  resonance  peaks  where  the  cross  section  values  are 
relatively  large.  In  these  vicinities,  the  relative  error  criterion  is  more  easily 
achieved  than  absolute  error  so  non-typical  relative  error  is  not  found.  However, 
tolerance  is  only  checked  at  the  midpoint  when  determining  suitability  of  an 
energy  mesh  point.  This  method  of  tolerance  criteria  allows  error  to  exceed 
tolerance  anywhere  besides  the  midpoint  that  was  evaluated.  As  expected,  the 
energy  location  of  non-typical  relative  error  is  where  the  cross  section  values  are 
near  zero.  In  this  case  relative  error  becomes  the  more  restrictive  criteria  and 
still  provides  no  accuracy  guarantees  at  points  other  than  the  midpoint. 

The  tolerances  specified  by  the  test  method  are  0.01  barns  for  absolute 
error  and  0.0001  for  relative  error.  Of  the  100,000  cross  section  values  evaluated 
99,980  have  an  absolute  error  less  than  the  absolute  tolerance  specified  and 
99,960  have  a  relative  error  less  than  the  relative  error  specified.  The  NJOY 
criterion  for  accuracy  is  not  documented.  So,  a  similar  quantified  statement 
cannot  be  made.  However,  box-whisker  plots  provide  a  means  to  compare  the 
extent  of  the  non-typical  error  of  NJOY  approximations  to  those  of  the  test 
method  approximations. 
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Box- Whisker  plots  depict  the  values  of  the  data  points  by  placing  them  on 


a  vertical  line  ranging  from  the  data  point  with  the  lowest  value  to  that  of 
highest  value.  The  data  points  found  between  the  first  and  third  quartile  of  the 
total  data  set  (middle  fifty  percent  of  the  data)  are  depicted  by  a  box  containing 
a  horizontal  line  that  is  the  median  [9].  A  whisker  at  the  value  where  data  points 
with  a  greater  value  are  considered  non-typical,  is  depicted  by  a  horizontal  line 
above  the  box  [9].  The  value  indicating  the  beginning  of  the  non-typical  data  is 
three-halves  the  value  marking  the  top  edge  of  the  box  (the  value  of  the  third 
quartile). 

Figure  15,  shows  the  box  whisker  plots  for  the  absolute  error  data  of  both 
NJOY  and  the  test  method.  For  ease  of  viewing  and  comparing  the  inter¬ 
quartile  ranges  and  whiskers  of  the  plots,  the  entire  range  of  the  data  is  not 
presented  in  Figure  15.  However,  quantified  values  for  the  graphical  depiction  of 
Figure  15,  as  well  as  other  relevant  statistical  values  not  displayed  are  provided 
in  Table  2.  From  Figure  15  it  is  seen  that  the  test  method  provides  a  more 
compact  data  set  than  NJOY.  The  whisker  indicating  the  location  of  outliers  is 
not  eliminated  but  is  closer  to  the  box  than  that  of  NJOY.  Information  provided 
in  Table  2  also  supports  the  test  method  as  more  consistent  approximation  by 
the  reduction  in  the  value  of  the  maximum  error  by  two  orders  of  magnitude  and 
more  importantly  the  reduction  in  non-typical  data  points  by  thirty-six  percent. 
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Figure  15.  Box-whisker  plot  of  absolute  error  for  NJOY  and  the  test  method 


NJOY 

TEST 

Mean 

0.00667 

0.00029 

Std.  Deviation 

0.08987 

0.00058 

Median 

0.00019 

0.00014 

1st  Quartile 

0.00007 

0.00004 

3rd  Quartile 

0.00055 

0.00036 

Inter-quartile  Range 

0.00048 

0.00032 

Non-typical  Value 

0.00128 

0.00082 

#  Non-typical  data  pts.  out  of  100,000 

12,093 

7684 

Max  Error  of  sample 

5.75882 

0.04107 

Table  2:  Statistical  data  of  absolute  error  (barns)  for  NJOY  and  the  test  method 
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Figure  16.  Box-whisker  plot  of  relative  error  of  NJOY  and  the  test  method 


NJOY 

TEST 

Mean 

0.00009 

0.00002 

Std.  Deviation 

0.00074 

0.00003 

Median 

0.00002 

0.00001 

1st  Quartile 

6.40851xl0"6 

4. 27699xl0"6 

3rd  Quartile 

0.00006 

0.00003 

Inter- quartile  Range 

0.00005 

0.00003 

Non-typical  Value 

0.00013 

0.00007 

#  Non-typical  data  pts.  out  of  100,000 

11,032 

4357 

Max  Error  of  sample 

0.12570 

0.00611 

Table  3:  Statistical  data  of  relative  error  for  NJOY  and  the  test  method 
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Figure  16,  shows  box  whisker  plots  for  the  relative  error  data  of  both 
NJOY  and  the  test  method.  Again,  for  easy  visual  comparison  of  the  inter¬ 
quartile  range  and  11011-typical  data  location,  the  entire  range  of  the  data  is  not 
presented,  but  quantified  values  are  provided  in  Table  3.  Figure  16,  shows  that 
the  test  method  provides  a  more  compact  data  set  than  NJOY  indicating  more 
consistent  results.  Information  provided  in  Table  3  supports  this  conclusion  by  a 
reduction  in  the  value  of  the  maximum  error  by  a  factor  of  two  and  a  reduction 
in  non-typical  data  points  by  sixty-one  percent. 

IV. 3:  Performance  Profiling  of  Linear  and  Cubic  Hermite  Approximation 

Comparisons  of  CPU  time  required  to  calculate  SLBW  resonance  cross 
sections,  piecewise  linear  approximations  and  piecewise  cubic  Hermite 
approximations  were  conducted  using  a  test  program  that  tracked  CPU  time  of 
the  respective  calculations.  The  energy  mesh  of  both  NJOY  and  the  test  method 
were  used  as  input  for  respective  calculations. 

Initial  comparisons  between  PAXK  linear  interpolation  and  piecewise  cubic 
Hermite  approximation  gave  the  counter-intuitive  result  that  piecewise  cubic 
Hermite  approximation  was  faster  than  linear.  This  result  was  found  true,  but  is 
attributed  to  the  logical  checks  required  by  PAXK  to  decide  which  type  of 
interpolation  to  execute,  be  it  two-dimensional,  logarithmic,  etc.  Because  only 
the  linear  interpolation  for  the  resonance  region  is  interesting,  these  logic  checks 
were  bypassed  for  the  comparison.  It  is  still  noteworthy  that  the  use  of  the  test 
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method  in  PAXK  applications  will  increase  performance  in  approximation  of  the 


SLBW  resonance  region  cross  sections. 

Table  4  provides  CPU  time  results  for  five  million  calculations  of  SLBW 
resonance  region  cross  sections,  piecewise  linear  approximations,  piecewise  cubic 
Hermite  approximations  as  well  as  the  sorting  time  required  to  locate  the 
appropriate  energy  mesh  location  required  to  calculate  the  respective 
approximation.  The  calculation  times  of  Table  4  were  measured  using  a  Pentium 
4,  3.2  gigahertz  processor  with  896  megabytes  of  available  random-access 
memory. 


CPU  TIME  IN  SECONDS 

SLBW  cross  section 

48.5156 

Linear  approximation 

0.12500 

Cubic  approximation 

0.20312 

Table  4:  CPU  time  requirements  for  SLBW  calculation  and  approximation 

The  time  required  to  calculate  piecewise  linear  approximation  to  SLBW 
cross  sections  is  three-tenths  of  a  percent  of  the  time  required  to  calculate  the 
same  SLBW  cross  sections.  Piecewise  cubic  Hermite  approximation  takes  four- 
tenths  of  a  percent  of  the  time  to  calculate  the  SLBW  cross  sections.  This  is  an 
increase  over  linear  approximation  of  one  tenth  of  a  percent.  However,  the  CPU 
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timing  results  indicate  that  both  methods  of  approximation  take  trivial  amounts 
of  CPU  time  when  compared  to  SLBW  calculations. 
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VI  Conclusion  and  Recommendations 


VI.  1:  Summary 

The  test  method  proved  adequate  in  meeting  the  goals  of  the  research  for 
the  SLBW  resonance  region.  The  test  method  reduced  the  number  of  energy 
mesh  points  by  a  factor  of  four  and  reduced  the  data  storage  requirement  of  the 
energy  mesh  by  a  factor  of  three.  Mean  accuracy  of  the  approximation  was 
improved  by  an  order  of  magnitude  for  absolute  error  and  a  factor  of  four  for 
relative  error.  Each  of  the  three  techniques  applied  by  the  test  method  is  simple 
and  easily  implemented  in  modularized,  modern  Fortran  90/95  code. 

The  use  of  piecewise  cubic  Her  mite  approximation  was  the  most 
significant  contributor  to  the  increase  in  accuracy  and  reduction  in  energy  mesh 
points  of  the  three  techniques  used  by  the  test  method.  The  adaptive  mesh 
calculator  was  the  weakest  performer  of  the  three  techniques  applied  by  the  test 
method  to  generate  the  resonance  region  cross  section  approximation.  The 
simple  bisection  adaptive  mesh  calculator  used  in  this  research  was  inferior  to  the 
methods  NJOY  uses  to  complete  this  task.  The  bisection  adaptive  mesh 
calculator  when  using  linear  interpolation  requires  an  increase  in  mesh  points  by 
a  factor  of  two  in  order  to  come  within  two  orders  of  magnitude  of  the  mean 
relative  error  of  NJOY. 

Use  of  the  cubic  Hermite  interpolation  increases  the  CPU  time 
requirements  slightly  over  use  of  linear  interpolation  in  the  general  sense. 
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However,  implementation  into  PAXK  increases  the  CPU  efficiency  of  the 
approximation. 

VI.2:  Recommendations  For  Future  Research 

The  reduction  in  mesh  points  of  the  test  method  from  that  of  the  NJOY 
mesh  proved  a  slight  advantage  in  total  approximation  CPU  time.  The 
advantage  was  the  result  of  the  search  routine  used  to  locate  the  energy  of 
interest  between  the  appropriate  two  points  in  the  energy  mesh.  Having  less 
mesh  points  to  check  improves  the  efficiency  of  the  search  algorithm.  This 
advantage  could  possibly  be  enhanced  to  the  point  where  piecewise  cubic  Hermite 
interpolation  to  could  exceed  the  CPU  efficiency  of  piecewise  linear 
approximation  if  a  more  sophisticated  adaptive  mesh  calculator  were  developed 
and  energy  mesh  points  were  further  reduced. 

A  similar  computational  advantage  to  Doppler  broadening  calculations  is 
also  expected  as  a  result  of  decreased  energy  mesh  points.  Doppler  broadening 
calculations  adjust  the  cross  section  values  from  0  K  to  temperatures  of  interest. 
Doppler  broadening  calculations  require  the  summation  of  integrals  over  adjacent 
points  in  the  energy  mesh  to  complete.  A  reduction  in  energy  mesh  points  is 
directly  proportional  to  a  reduction  in  integrals  required  for  Doppler  broadening. 
These  integrals  are  non-trivial.  So,  the  reduction  in  the  number  of  integrals 
required  to  perform  Doppler  broadening  has  the  probability  of  significantly 
decreasing  the  CPU  calculation  time  for  Doppler  broadening  calculations  if 
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piecewise  cubic  Hermite  approximation  is  used  instead  of  piecewise  linear 
approximation. 
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Appendix  A:  Nuclear  Cross  Sections 


Nuclear  cross  sections  represent  the  probability  of  particle  interactions 
with  the  nucleus.  The  microscopic  cross  section  has  units  of  area,  specifically 
barns,  where  lbarn  =  10-28 meters2 .  The  cross  section  includes  more  than  the 
physical  area  of  the  target  nucleus.  Quantum  mechanics  influences  the 
probability  of  particle  interaction  and  changes  the  effective  area  of  the  nucleus. 
Cross  sections  often  depend  greatly  on  the  energy  characteristics  of  bombarding 
particles  as  well  the  quantum  characteristics  of  the  nucleus  [10]. 

There  are  many  neutron-nucleus  interaction  types.  These  include  elastic 
scattering  (which  is  of  interest  in  the  main  text),  inelastic  scattering,  neutron 
capture  and  fission.  The  first  two  types  are  combined  to  form  the  scattering 
cross  section  while  the  latter  two  types  combine  to  form  the  absorption  cross 
section.  The  absorption  and  scattering  cross  section  combine  to  form  the  total 
cross  section  [3]. 

Cross  sections  can  also  represent  the  probability  of  interaction  per  unit 
length  of  travel  in  material  as  opposed  to  interaction  with  a  single  nucleus.  This 
probability  is  the  macroscopic  cross  section.  The  macroscopic  cross  section  is 
related  to  the  microscopic  cross  section  by 

-  W{«2S“  M— [11].  0) 

meters  atom  barn 
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The  mean  free  path  is  the  average  distance  traveled  in  a  material  before 


interaction.  The  inverse  of  the  macroscopic  cross  section,  has  units  of  length  and 
gives  the  value  of  the  mean  free  path  of  the  neutron  (or  other  bombarding 
particle  such  as  a  photon)  [3]. 
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Appendix  B:  ENDF  Cross  Section  Information 


Evaluated  Nuclear  Data  Files  version  B-VI,  ENDF /B-VI,  are  libraries  that 
contain  computer  readable  files  that  describe  nuclear  reaction  cross  sections,  the 
distributions  in  energy  and  angle  of  reaction  products  among  other  properties 
[12].  The  ENDF  libraries  contain  a  collection  of  documented  data  evaluations 
that  can  be  used  as  the  primary  input  to  nuclear  data  processing  codes.  ENDF 
files  are  specific  to  a  material.  A  material  can  be  single  nuclide  (isotope),  a 
natural  element  containing  several  isotopes,  or  a  mixture  of  several  elements 
(compounds,  alloy  or  molecule)  [7].  The  ENDF  that  is  specific  to  the  single 
isotope  Pu238  was  used  as  the  source  of  parameters  for  calculations  in  the  main 
text. 

The  ENDF  for  a  material  can  contain  up  to  twenty-six  sub  files  each 
containing  specified  information  on  various  nuclear  and  cross  section 
characteristic.  The  Pu238  ENDF  only  contains  six  of  the  possible  twenty-six 
files.  Of  these  six,  only  File  2  is  required  for  the  calculation  of  resonance 
contributions  to  the  cross  section.  File  2  provides  all  required  parameters 
necessary  to  evaluate  resonance  contributions  to  the  cross  section.  These 
parameters  include  all  those  required  for  calculation  of  the  SLBW  formula  of 
Appendix  C  as  well  as  additional  identifiers  such  as  the  left  and  right  energy 
boundaries  of  the  resonance  region.  For  the  case  of  SLBW  calculations  these 
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contributions  include  the  elastic  scattering,  radiative  capture  and  fission.  Any 
other  contributing  reaction  to  the  total  cross  section  in  the  resonance  region  is 
provided  in  File  3  and  added  to  the  calculations  using  File  2  parameters  for  a 
given  energy.  File  3  contains  information  in  the  form  of  energy  cross  section 
pairs  with  parameters  that  specify  the  interpolation  method  used  to  calculate  the 
desired  cross  section  for  a  given  energy.  For  the  case  Pu238  ENDF  provides  no 
File  3  contributions  to  the  cross  section  in  the  SLBW  resonance  region. 
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Appendix  C:  The  SLBW  Elastic  Scattering  Formula  and  First 

Derivative 

The  formula  for  SLBW  elastic  scattering  cross  section  requires 
examination  of  the  documentation  from  ENDF-B/VI  for  basic  form.  In  order  to 
obtain  a  form  useful  in  computation,  examination  of  the  NJOY  source  code  is 
also  required.  The  presentation  of  the  elastic  scattering  cross  section 
formalization  that  follows  is  one  not  directly  found  in  ENDF  documentation  but 
implemented  and  validated  by  both  NJOY  and  PAXK.  The  derivation  of  the 
Derivative  was  a  critical  part  of  this  research  for  which  no  other  reference  was 
used. 

C.l:  SLBW  Resonance  Region  Formalism 
The  elastic  scattering  cross  section  is, 

nls— 1 

(E)  =  5]  <x,„  (E).  (4) 

1= 0 

Where,  l  is  the  neutron  angular  momentum  number  and  nls  is  an  integer  in 
ENDF  that  equates  to  the  number  of  l  values  considered  for  a  particular 
resonance  region.  Currently,  in  ENDF  files,  The  maximum  nls  value  is  four [7]. 

The  first  expansion  of  the  elastic  scattering  cross  section  is  given  in 
Equation  Five  on  the  following  page, 
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aell(E)  =  (2l  +  l)—^m2(4>i(E))  + 


Jmax  NR, 


h  sr'  _ 

fc2  (E)  ,6l 9jh[{E  -  E'r  (E))2  +  y4T2r  (E) 


Kr  (E) cos  (2^  (F))  -2rn>r  (£7)r,ir  sin2  (F)) 
2(^-E;(E))rriirsin(20i(E)))). 


Where, 


h„  =  l  +  I  H - 


2J  +  1 
2(2/  + 1) ' 


The  parameter  I  is  the  target  spin  number  and  NR  is  the  number 


of  resolved  resonances  for  a  given  l  -state. 


C.2:  Sub-Functions  and  First  Derivatives 


The  neutron  wave  number  in  the  center  of  mass  system  (“A”,  is  the  ratio 
of  the  mass  of  a  particular  isotope  to  that  of  the  neutron): 


k(E)  =  0.00219680323015523  * - Vf 

A  +  l 


/'(F) =  0.00219680323015523 


A  1 
A  +  l  2>/F  ' 


The  quantity  pis, 


p(F)  =  k{E)AP  -►  p  (E)  =  k\E)AP. 


Where,  AP  is  the  scattering  radius. 
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The  quantity  p  if  NAPS  =  0  (NAPS  is  a  parameter  that  controls  use  of 


the  channel  radius  a  and  the  scattering  radius  AP)  is, 

p{E)  =  k(E)a  ->•  p\E)  =  k'(E)a 
p  =  k{E)AP  -►  p(E)  =  k\E)AP 

The  quantity  if  NAPS  =  1  is, 

p(E)  =  p{E)  =  k{e)AP  ->•  p(E)  =  p{E)  =  k\e)AP. 

Where,  AP  is  provided  in  the  material’s  ENDF  and, 

i 

a  =  0.123(A  *  1.00866491578)3  +.08 


(10) 


(11) 


(12) 


The  (negative  of  a)  hard-sphere  phase  shift  is, 


AW  =  p(E)  ->  A'  =  p’ W 

<p2(E)  —  p(E)  —  arctan  (p(l?))  — +  A  —  p1  (E)  — 


P  (E) 


1  +  P  W 


AW  —  p(E)  —  arctan 
<j>2{E)  —  p(E)  —  arctan 

A'W  =  />'W* 


3  P(E) 


3-/5  W 

p(E)(lb-p(E)) 


A  =  AW 


P‘(E) 


9  +  3 p2  (E)  +  p‘  (E) 


15  -  6p2  (E) 


(13) 


p‘  (E)  -  60£  (E)  +  1392/5*  (E)  -  13500/53  ( E )  +  49050^  ( E )  +  900p(E)  -  3150 
p‘  (E)  -  60 p‘  (E)50445p2  (E)  +  1386p4  (E)  -  13500p3  (E)  +  50445 j?  (E)  +  225  ' 


The  penetration  factor  is  given  by  Equation  Fourteen  on  the  following 


page, 
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P„(E)  =  P(E)  P'  =  p'(E) 


P,(E) 

P,(E) 

P,(E) 

P,(E) 


PM  .  3p](-E)  2+£)  ' 

1 +  p2(E)  1  +  +  p2(£)  (1  +  p’(E)f’ 

P‘(E)  ,  p,(E){pt(E)+9p2(E)  +  45)  (14) 

9  +  3 p2(E)  +  p'(E)  2  H  (p4(£)  +  3p2(£)  +  9)2 

_ P+) _ , 

225  +  45  f>\E)  +  6  p\E)  +  p\E) 

f/‘(E)(p‘(E)  +  18  p\E)  +  225p2(g)  +  1575) 

(p6(E)  +  6p4(E)  +  45p’(E)  +  225)2 


The  neutron  width  is, 


r  (F)  P^E)GNt  .  P;(E)GNr 

"A  ’  j=;(|£’i?D  ' 

Where,  GNr  is  the  neutron  width  at  Er  provided  by  ENDF. 


(15) 


The  total  width  is, 


r,.(£)  =  r  (£)  +  r  ->  r.)  =  r  ’(£). 


(16) 


Where,  Yxr  is  a  sum  of  the  ENDF  provided  parameters  GG  (the  radiation 
width)  and  GF  (the  fission  width).  The  shift  factor  is, 


S0  =  0 


Si(E) 


-1 


9  -*S[(E)  =  p\E)  2pi-E)  9, 

1  +  P\E)  (l  +  P\E)f 


18  +  3  p2{E) 


S  iE]  ~  ^  v+z  :  S,(E)  6^)(9  +  12p2(F)  +  p4(F)) 

^2\E)  n  „  2/tp\  ,  4/n\  ^2\E)  P  ,  „  ,  \2  ’  ^  ' 


Ss(E)  =  ~ 


9  +  3  p\E)  +  p\E) 

675  +  90p2(F)  +  6p4(F) 

225  +  45p2(F)  +  6p\E)  +  p6(E) 


9  +  3  p\E)  +  p\E) 


6p(F)(3375  +  1800p2(E)  +  765p4(F)  +  60p6(F)  +  2  p8(E) 

bn(E)  —  p  (E)  - 

225  +  45p2  (E)  +  6p 4  (E)  +  p\E) 
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The  primed  resonance  energy,  if  /  equals  zero  is, 


E'r  =  Er  -t  \E'r  1  =  0. 


(18) 


If  /  equals  one  the  primed  resonance  energy  is, 


±jr  ±jr  ~r  _ ,i  „  k  -1-  n,r\\-^r\ 


(19) 


m\ET 


C.3:  The  First  derivative  of  <jel 

The  derivative  of  the  SLBW  elastic  scattering  nuclear  cross  section  is, 


nls— 1 

<1  (E)  =  £  ( E )  (20) 

/=0 

In  order  to  show  the  form  of  <j'eI  in  its  first  expansion  it  is  useful  to  accumulate 
the  terms  of  equation  (6)  to  produce  the  following  equation, 


NR, 

ael,l  (E)  =  f(E)  +  G(E)  £  9.7  £'!(£)• 

r=l 


(21) 


“<j'el”  then  becomes, 

Jmax  Jmax  NR, 

<,  ( E )  =  f'(E)  +  G\E)  £  <,.,£*(£)  +  <?(£)  £  s,£fcW-  (22) 

J—Jr  nin  1  J  — J,ni„  f— 1 


Where, 


/(£) 


8(2/  +  l)7r  cos  (0;  (T1))  sin  (A1))  (/>,'  (E) 


8(2/  +  l)7r  sin2  (A1))  A:'  (£/) 


.(23) 


And, 


50 


2nlz  (E) 

W 


(24) 


G'  (E)  = 


And, 

ft  =  $  QE)  r„,  ( -E)  (4  cos  (20,  (£))(£-  E’r  (£Q)  -  2  cos  (2^  (£Q)  r„,  (g)  _ 

(£--B,'(£))2  +  ^r;(£) 

(g)r»,,-(£)(4  costal  (-E))sin(,)>,  (-E))r,,)  | 

(B-B,'(B))2  +  i^r).(B) 

2r,,sin(2^(g))r,r(£)(l-[g,'],(g)) 

(£-B,'(iJ))s+%r)(B) 

UgH sin(2^,,  (£))  (1  -  g;  (£Q)  +  2  cos(2fl  (£))r„,(g)  -  2r, ,  sin2  (0,  (£))) 

(£-£,'(£))2+^r)(B) 

(r„,(£)(r„,  ( E)coS{2<p,  (E))  -  2Txr  sin2  (£))  +  2{E  -  E'r  (E)) sm(2^  (£))))  * 
(2(1  -  E'r  (E))(  1  -  E'r  (E))  +  irr(£)rr'(£)) 

((£-£,' TO)2 +  Mr;W  ’ 


(25) 
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Appendix  D:  The  Shifted  and  Scaled  Cubic  Hermite  Polynomial 

The  cubic  Hermite  interpolating  polynomial  is  of  the  form, 

f(x)  =  c0  +  x{cx  +  x{c2  +  c3x)).  (26) 

Equation  24,  has  un-scaled  or  un-shifted  coefficients: 

f{b)a2{a  —  36)  +  ( 6(/'(6)a2(6  —  a)  +  6(/(a)(3a  —  b)  +  f'(a)a(b  —  a) 

c°  =  7  7\3 

(a-b) 

_  /'(6)q3  +  o(6/(6)  -  6 f(a)  +  (2 f\a)  +  f(b))a)b  -  (/'(a)  +  2 f\b))ab2  -  f(a)bs 
1  (a  -  bf 

3/(a)(a  +  b)  —  3/(6)  (o  +  b)  —  (a  —  &)(/  (o)o  +  2/  (6)a  +  2/  (o)6  +  /  (6)6) 

C2  _  (a  -  6)3 

_  2/(6)  -  2/(a)  +  (/'(a)  +  /(6))(a  -  6) 

3  (a  -  bf 

Where,  a  and  b  are  the  endpoints  of  the  cubic  Hermite  polynomial. 

Unfortunately  these  coefficients,  although  correct  will  not  work  when  using 
the  adaptive  mesh  calculator.  This  is  especially  apparent  in  the  vicinity  of  the 
resonance  peaks  where  mesh  points  get  extremely  close  together.  The  cubed  term 
in  the  denominator  becomes  very  small  as  a  and  b  get  very  close  together.  This 
results  in  catastrophic  cancellation  in  the  denominator  and  overflow  of  the 
coefficients. 

By  scaling  the  function  to  where  a  and  b  always  have  a  fixed  value,  the 
catastrophic  cancellation  and  overflow  is  prevented  by  fixing  the  distance 
between  the  endpoints  of  the  cubic  function  to  reasonable  value.  The  change  of 
variable  that  leads  to  this  desired  effect  is  easiest  to  accomplish  when  the  scaled 
function  is  also  shifted  to  be  symmetric  about  zero.  The  derivation  of  the  scaled 
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Hermite  Cubic  and  coefficients  are  as  follows:  Set  the  function  u(x)  such  that 

when  x  =  a  u  =  —  —  and  when  a;  =  b  =$>  u  =  — .  Under  these  conditions 

2  2 


u(x)  becomes, 


x  — 


b  +  a 


u{x)  = 


b  —  a 


(28) 


Now, 


g(u(x))  =  k0  +  u +  u(k2  +  k3u)).  (29) 

Setting  g(u(x ))  equal  to  f(x)  of  equation  (30)  and  solving  for  the  coefficients  of 


g(u(x))  provides, 


+  m  ,  (6-a)(/'(a)-/'(6)) 


+ 


K  = 


k,  = 


3  (f(b)~f(a))  {b-a){f'{a)  +  f\b)) 

2  4 

(f\b)  -  f'(a))(b  -  a) 


(31) 


^  =  2(/(a)  -  /(6))  +  (6  -  a)(/'(a)  +  /'(&)). 


These  scaled  coefficients  were  used  to  produce  the  piecewise  cubic  Hermite 
polynomial  cross  section  representation  presented  in  the  main  text. 

It  is  important  to  note  the  potential  loss  of  accuracy  when  executing  the 
change  of  variable  from  v  to  u(x)  in  computing  applications.  The  accuracy  of 
u(x)  is  reduced  from  that  of  x  in  any  region  other  than  where  the  true  value  of 
a  and  b  is  negative  one-half  and  one-half  respectively  (a  condition  that  is  never 
true  for  the  application  of  the  scaled  cubic  used  in  the  main  text).  This  loss  of 
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accuracy  has  potential  to  be  significant  in  the  regions  where  the  un-scaled 
coefficients  experience  catastrophic  cancellation  and  overflow.  When  the  range 
between  a  and  b  is  small  the  value  of  w(x)will  most  certainly  be  much  larger 
than  the  value  of  x  reducing  the  number  of  available  digits  to  the  right  of  the 
decimal  place.  This  loss  of  accuracy  will  propagate  through  the  evaluation 
of g(u{x)). 
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Appendix  E:  Algorithms  for  SLBW  Resonance  Region  Cross  Section 

Approximation 

E.l:  Calculation  of  the  Initial  Energy  Mesh 

Input:  Left  and  right  boundaries  of  resonance  region  and  Er  values. 

Output:  Initial  energy  mesh  of  size  2 rm£LX  +2  containing  left  and  right 
boundaries  of  the  resonance  region  and  local  minima  and  maxima  of  the 
resonance  peaks. 

For  all  Er ,  r  =  1 . rmax , 

a  =  Er  —  A  E 

b  —  Er  +  A  E. 

While  sign  a' (a)  is  negative  or  sign  a\b )  is  positive, 

a  =  a  —  A  E 
b  =  b  +  AE. 

When  sign  cr' (a)  is  positive  and  sign  a \b)  is  negative, 

Bisection  root  solve  for  E  where  cr'(E)  =  0  ,  between  a  and  b  . 

Then, 

77  _  77 

J~Jr, local  a  max 

b  =  a 

a  =  a  —  2A  E. 

While  sign  cr'(a)  is  positive, 
a  =  a  —  A  E  . 

When  sign  a  '(a)  is  negative  and  sign  a\b)  is  positive, 

Bisection  root  solve  forE  where,  a1  (E)  =  0  is  between  a  and  b  . 
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Then, 


E  =  E 


r, local  a  min 


E.2:  Adaptive  Energy  Mesh  Calculator 


Input:  Sorted  initial  energy  mesh  oiEn,  where,  n  =  1.2,  3 . initial  max-1, 

initial  max. 

Output:  Energy  mesh  of  required  size  to  meet  accuracy  requirements. 


Starting  with  n  =  2  , 

midpoint=  — — 

2 


'  approximated 


(midpoint)  is  within  tolerance  of  <x( midpoint) ,  then  add  “midpoint’ 


to  the  final  energy  mesh. 


'  approximated 


(midpoint)  is  not  within  tolerance  of  er(midpoint) ,  then 


midpoint= 


midpoint  —  En 


’  approximated 


(midpoint)  is  within  tolerance  of  a(midpoint) ,  then  add  “midpoint’ 


to  the  final  energy  mesh  and, 

midpoint=— — m^P°^n^  an(]  repeat.  above  process. 

2 

If  midpoint=  — — niidpoint  ^  j-j.ie  f}naj  energy  mesh  then, 

2 

E  ,  —  E 

midpoint—  - —  and  repeat  above  process. 


When,  midpoint— 


—  midpoint  . 


is  added  to  the  final  energy  mesh  exit. 
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E.3:  Piecewise  Shifted  and  Scaled  Cubic  Hermite  Approximation 

Input:  Left  boundary  (a),  right  boundary  ( b ),  energy  ( E )  between  a  and  6, 
a  (a),  cr(6),  cr' (a)  and  cr'(b). 

Output:  Cubic  Hermite  approximation  of  cr{E) .  Uipprox(A) 

Calculate  the  piecewise  cubic  Hermite  coefficients 

k  _(cr(a)  +  CT(b))  |  ( b-a)(a'(a)-a'(b )) 

0  2  8 
j  _3 (<j(b)-<j(a))  (b  -  a)(<j' (a)  +  <j\b)) 
l~  2  4 

k  _  (a'(b)-a'(a))(b-a) 

2  2 

£3  =  2(cr(a)  -  cr(6))  +  (b  -  a)(<j'(a)  +  <j’(b )). 

Change  energy  variable  to  the  scaled  polynomial, 

E-^ 

u= - 2-. 

b-a 

Calculate  the  cubic  Hermite  approximation, 

^approx  (E)  =  k0  + 11  (k\  +  u(k2  +  k2u))- 
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